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We have combined ab initio quantum chemistry calculations with a rate-equation formalism to analyze elec- 
troluminescence spectra in single-molecule junctions, measured recently by several groups in scanning tunnel- 
ing microscope setups. In our method, the entire vibrational spectrum is taken into account. Our method leads 
to good quantitative agreement with both the spectroscopic features of the measurements and their current and 
voltage dependence. Moreover, our method is able to explain several previously unexplained features. We show 
that in general, the quantum yield is expected to be suppressed at high bias, as is observed in one of the mea- 
surements. Additionally, we comment on the influence of the vibrational relaxation times on several features of 
the spectrum. 

PACS numbers: 78.60.Fi, 31.15.A-, 73.63.-b, 85.65. +h 




FIG. 1. Molecular structure of the isolated porphyrin derivatives 
used in the calculations: (a) ZnEtioI, measured by Qiu et ah (Ref. 51) 
and (b) H2TBPP, measured by Dong et al (Ref. 52). 



I. INTRODUCTION 

Scanning tunneling microscopy (STM) has proven to be 
an invaluable tool for studying single molecules and their 
potential applications in nanoscale devices. In addition 
to measuring conductance properties, the electrolumines- 
cence spectrum, induced by the tunneling current can be 
studied."^"^ Electroluminescence has been experimentally ob- 
served on clean metal surf aces, films,^^'^^ nanoparticles,^^ 
and single molecules.^ ^"^"^ Additionally, electroluminescence 
has been observed in lithographically defined three-terminal 
CdSe single-nanocrystal junctions. In the case of single 
molecules, the simultaneous measurement of the electrical 
and optical behavior has the potential to greatly enhance our 
understanding of nanoscale junctions, and, through detailed 
analysis of the vibrational spectrum, to provide valuable in- 
sight into the conformational structure of single molecules in 
a junction.^^ 

We have developed a method based on a combination of 
ab initio density functional theory (DFT) calculations with a 
rate-equation formalism to analyze electroluminescence mea- 
surements on weakly coupled single-molecule junctions. In 
this paper we apply our method to two STM measurements on 
single porphyrin derivatives: Zn(II)-etioporphyrin I [ZnEtioI, 
see Fig. 1(a)] measured by Qiu et al.,^^ and m^^(9-tetrakis(3,5- 
diterbutylphenyl)porphyrin [H2TBPP, see Fig. 1(b)] mea- 
sured by Dong et al^^ It is known that the relative elec- 



tronic coupling strengths to the source (S) (STM tip) and 
drain (D) (substrate) electrodes greatly aff'ect molecular 
electroluminescence,^^ and in order to suppress quenching 
due to the substrate, it is necessary to employ a spacer layer. 
Qiu et al^^ have used a thin AI2O3 film on the NiAl(llO) 
substrate to act as a spacer layer. Although the electronic 
coupling to the leads is still highly asymmetric, electrolumi- 
nescence is sufficiently enhanced to be observable. Dong et 
al.,^^ on the other hand, have deposited several monolayers 
of H2TBPP, resulting in a nearly symmetric coupling of the 
molecule to the leads. Since these measurements investigate 
similar molecules in diff'erent regimes (asymmetric and sym- 
metric coupling), they provide a good test case to understand 
the physics of single-molecule electroluminescence. 

It is known that plasmons in the substrate and the STM tip 
can mediate photoemission, and in the case of surface en- 
hanced Raman scattering this can even be a powerful spec- 
troscopic tool.^^ However, in order to observe electrolumi- 
nescence originating solely from a single molecule, the sig- 
nal from plasmons should be suppressed as much as possible. 
In the experiments under discussion, several measures have 
been taken to this end, including using atomically flat sub- 
strates and tungsten STM tips. Also, in the measurements of 
Qiu et al,^^ applying a series of high- voltage pulses between 
the tip and the substrate seems to make the plasmon spectrum 
smoother. Although an enhancement of the photoemission 
rate due to the STM tip and the substrate probably still ex- 
ists, the experiments show that this enhancement is generally 
rather structureless and does not afl'ect the shape of the spec- 
tra. Moreover, control measurements in both groups show that 
the electroluminescence intensity of the bare surface is much 
weaker than that of the molecules.^ ^'^^ We will therefore ig- 
nore the efl'ect of plasmons from now on. 

The measurement of Qiu et aO^ has been previously 
analyzed by Buker and Kirczenow using the Lippmann- 
Schwinger Green's function scattering technique. ^^'^^ In this 
paper we employ the rate-equation formalism, which allows 
us to take the entire vibrational spectrum of the molecule into 
account. 
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FIG. 2. Schematic picture of sequential tunneling through a molec- 
ular junction. Fs and Fd are the electronic couplings to the source 
and drain electrodes, respectively, and a is the electrostatic coupling 
to the bias voltage. These quantities can, in principle, be different for 
different orbitals. The thick solid lines represent the HOMO of the 1^ 
and neutral charge state, and the dashed line represents the LUMO 
of the neutral charge state. The thin lines are vibrational excitations. 
The small vertical arrow at the top indicates vibrational relaxation, 
and the wavy line indicates the emission of a photon. 




II. METHOD 

The magnitude of the electronic coupling constants ob- 
tained from the total current and the JZ/JV-curves in the mea- 
surements (discussed in Sec. Ill A) is indicative of Coulomb 
blockade [see Fig. 2 in Ref. 51, and Fig. 5(a)], suggesting that 
the molecule-electrode coupling is weak (Fs, Fd, k^T ^ AE), 
and the electron addition energies (AE) only allow the tun- 
neling of one electron at a time (sequential tunneling). A 
schematic picture of this process is shown in Fig. 2.^^ In 
this picture, current can flow as soon as the highest-occupied 
molecular orbital (HOMO) of the neutral charge state enters 
the bias window. When also the lowest-unoccupied molecu- 
lar orbital (LUMO) becomes available, it is possible for the 
molecule to be in an electronically excited state when an elec- 
tron tunnels onto the LUMO instead of the HOMO. If the cou- 
pling to the drain (Fd) is weak enough, and the electron stays 
on the molecule for some time, the excited state can decay 
to the ground state and emit a photon. The efficiency of this 
process is given by the luminescence quantum yield, which 
is defined as the number of emitted photons per transmitted 
electron. 

In the weak-coupling regime, the current- voltage character- 
istics and the electroluminescence spectra can be calculated 
with the rate-equation formalism.^^"^^ In this formalism, a fi- 
nite number of molecular states is taken into account, and all 
processes are described in terms of transitions between these 



FIG. 3. The one-, two-, and three-electron states of a single 
molecule, with an even number of electrons in the neutral state and a 
non-degenerate HOMO and LUMO. Vibrational excitations are not 
shown. The labels are of the format ^"^^^Qm,, where Q is the charge 
of the molecule, S the total spin, and the eigenvalue of S^. Elec- 
tronically excited states are denoted by \ Transitions from the 1^ to 
neutral state are shown in red, transitions from the neutral to 1~ in 
blue, and the photoemission transition is shown in green. Ef^^~^ is 
the electron addition energy, E^^ the photo-emission energy, and £ex 
the exchange energy between the singlet and triplet state of 0'. Note 
that the energy differences are not to scale. 



states at a certain rate. The central quantity in this formalism 
is the vector of occupation probabilities Here n is taken 
to be a general quantum number describing charge, spin, and 
electronic and vibrational excitations. The time evolution of 
the occupation probabilities is governed by the master equa- 
tion. 



^ = y(P^W^ -PW 



(1) 



n' +n 



where y^n'^n and ^n^n' are the transition rate constants be- 
tween the states n and n' (given by Fermi's golden rule). This 
equation can be written in matrix- vector form. 



at 



(2) 
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where W is the rate-constant matrix with elements 



if / = j, 
otherwise. 



(3) 



The stationary states, with ^ = 0, correspond to the null 
space of W, with the condition that all elements of P are 
non-negative, and 2]„ = 1. 

The states of a single molecule with an even number of 
electrons in the neutral state, and a non-degenerate HOMO 
and LUMO are shown in Fig. 3. Ignoring vibrational excita- 
tions for the moment, the spin multiplicities of three charge 
states and one excited state result in a total of nine states, and 
therefore a 9 x 9 rate-constant matrix. However, assuming 
states with the same total spin have equal electronic (Fs and 
Fd) and electrostatic (a) (Ref. 72) couplings to the leads, and 
therefore equal transition rates, and ^1^^ can be com- 
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bined into ^1+, ^0^ , ^0' , and ^0' , into ^0', and ^IT and ^l" , 



^1' 



-1 



into ^1", leaving only five states. Since we are now describ- 
ing the system in terms of combined states, the transition rates 
diff'er from the original ones in that they have to be multiplied 
by a prefactor depending on the multiplicity of the states in- 
volved. Calling Mn = 2S n I the multiplicity of state n, the 
prefactor for the transition rate constant Wn^n' in this system 
is given by 



1 1 otherwise. 



(4) 



This also carries over into the case where vibrational ex- 
citations are taken into account, as long as they are spin- 
independent. 

The transition rates themselves depend on the type of tran- 
sition: tunneling, photo-emission, or vibrational relaxation 
(indicated by the arrows in Fig. 2). Tunneling of electrons 
to and from the source and drain is responsible for transi- 
tions between states with a diff'erent number of electrons on 
the molecule. Here we only consider single-electron tun- 
neling events. In the case of charging, where state n' has 
one more electron than state n, the rate constant is given by 
Wn^,r = WS_,„, + W^^„„ where 



= C F f 



(5) 



and correspondingly for W^_^^,, where Fnn' are the Franck- 
Condon factors (the overlap between the nuclear wave func- 
tions) (Refs. 71 and 73) and / is the Fermi function. Similarly, 
in the case of discharging, the transition rate constant is given 
by Wn'^n = W^.^^ + where 



and correspondingly for W^^„. Photo-emission is possible 
from the singlet excited state ^Oq to the ground state ^Oq (fluo- 
rescence). Transitions from the triplet states ^Op ^Oq, and ^0^^ 



(phosphorescence) are spin forbidden, and can only occur in 
the presence of spin-orbit coupling. Since phosphorescence 
occurs on a much slower time scale than fluorescence — and 
especially considering the fact that excited electrons can tun- 
nel ofl' the molecule in a junction — we will only take photo- 
emission from the singlet state into account. The radiative 
transition rate constant from the excited state n' to the ground 
state n is given by^"^ 



W - F / 



(7) 



where oj is the frequency of the emitted light (determined ap- 
proximately by the HOMO-LUMO gap), and yu is the transi- 
tion dipole moment. Note that we are only considering ra- 
diative transitions between the excited state and the ground 
state. Since the time scale of both the radiative and non- 
radiative transitions (microsecond to nanosecond) is typically 
much slower than the charging and discharging time scale (ap- 
proximately picosecond), these transitions will have a negli- 
gible efl'ect on the occupation probabilities. Only the radia- 
tive transitions are directly observable, while the non-radiative 
transitions will usually have an imperceptible efl'ect on the to- 
tal current. We will therefore only include the former in the 
rate equations and ignore the latter. Vibrational relaxation is 
taken into account with a single relaxation time for all vibra- 
tionally excited states. Separating the vibrational quantum 
number v from the other quantum numbers n, the transition 
rate constant for this relaxation is given by 



1 peq 

T 



where r is the relaxation rate and 



(8) 



(9) 



is the equilibrium occupation of the vibrational excitations ac- 
cording to the Boltzmann distribution. Note that a single re- 
laxation time for all vibrational excitations is only a good ap- 
proximation when relaxation is either much faster or much 
slower than all other rates. We will come back to this in 
Sec. IIIB. 

For the molecules under investigation, we have calculated 
the equilibrium geometry and the vibrational modes of the iso- 
lated system using DFT.^^ Although DFT is well suited for 
ground- state calculations, taking excited states into account 
can be problematic. We have therefore approximated the ex- 
cited state by forcing the molecule into a high-spin (5" = 1) 
configuration. Since the equilibrium geometry and the nor- 
mal modes only depend on the occupation of the molecu- 
lar orbitals through the charge density, and not on the total 
electron spin, this is not expected to have a significant ef- 
fect on the vibrational spectrum. The Franck-Condon factors 
have been calculated from the equilibrium geometries and the 
normal modes by using the method of Ruhoff' and Ratner,^^ 
while the transition dipole moments have been calculated with 
time-dependent DFT. This leaves the electronic and electro- 
static couplings, and the position of the Fermi level within 
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D 



(b) 
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FIG. 4. (a) Sequential tunneling with one charge state (0) and one 
excited state (0') in the bias window, (b) Sequential tunneling with 
two charge states (0 and 1~) and one excited state (0') in the bias 
window. The horizontal lines represent transitions between many- 
electron states. 



the HOMO-LUMO gap as fit-parameters. Note that these pa- 
rameters depend on the unknown contact geometry and vary 
from sample to sample in the measurements. Since extract- 
ing the HOMO-LUMO gap from DFT calculations is deli- 
cate, this value is fitted to the measurements as well. These 
fit-parameters only have a small influence on the vibrational 
features of the spectrum, which is the focus of our research. 

Before applying our method to the ZnEtioI and H2TBPP 
porphyrin derivatives, let us first consider simple situations 
which can be solved analytically. We ignore vibrational ex- 
citations for the time being. In Fig. 4(a), the and 0' states 
are in the bias window and available for transport, while the 
1" state remains unoccupied. In Fig. 4(b), also the 1" state is 
available. In both cases we assume the electronic and electro- 
static couplings to the leads to be equal for all states, and the 
temperature to be low enough for the Fermi functions to be 
either 1 or 0. Since both ZnEtioI and H2TBPP have a (nearly) 
degenerate LUMO, we double the multiplicities of the states 
involving the LUMO, resulting in the following five states: 
^1^, ^0, ^0', ^0', and ^1". As a final approximation, we assume 
the photo-emission rate to be much smaller than the charging 
and discharging rates. It will therefore have a negligible eff'ect 
on the occupation probability of the singlet excited state and 
can be omitted from the master equation. Since the rates for 
the singlet and triplet excited states are now equal, they can be 
combined into one excited state ^0', provided they are both in 
the bias window, leaving only four states to consider: ^1^, ^0, 



^O^andn-. 

In the case of Fig. 4(a), where the "^1" 
the rate-constant matrix is: 



state can be ignored. 



srs 2rD Fd 

Ts -2Tu 

1 4rs -FdJ 



(10) 



resulting in the stationary occupation probability 

1 



P = 



9rs + 2rD 



(2ru] 
Ursj 



(11) 



The current is given by 

n n 9ls + zFd 

and the electroluminescence intensity by 

1 2rsrE 



4 n 



7i9rs + 2rD 



(12) 



(13) 



where Te = tiW^, and the factor | is due to the fact that only 
the singlet, i.e., one quarter of the states in ^0', can emit a pho- 
ton. This equation shows that electroluminescence can indeed 
be quenched when the coupling to the drain is too large. The 
quantum yield is then 



Q 



I 



In the case of Fig. 4(b), the rate-constant matrix is: 



W 



(-5Ts 

Fs -4rs - 2Tj, 
4Ts 
4rs 



Ld 
Fd 

-Fs-Fd 2Fd 

Ts -3Fd 



(14) 



(15) 



resulting in the stationary occupation probability 



P = 



1 



(4Fs+Fd)(Fs+2Fd) 



2F^ 
TsFd 
SFsFd 



(16) 



The current is now given by 



/= rrs(5Pi. +4Po + /^oO 
n 

^ e 2FsFd(6Fs + 5Fd) 
;i(4Fs+FD)(rs+2FD)' 
and the electroluminescence intensity by 



L = 



iLe 

4 n 



2FsFdFe 



(4Fs+Fd)(Fs+2Fd)' 



resulting in a quantum yield of 

Q = 



Le 



6Fs + 5Fd 



(17) 



(18) 



(19) 



Comparing Eqs. 14 and 19 shows that the availability of a 
second charge state in the bias window changes the quantum 
yield by a factor of f2f^ff^' ^1 always decreases, with the 
magnitude of the change being determined by the asymmetry 
in the coupling. Although the expression for the change in the 
quantum yield depends on the particular multiplicities of the 
states involved, this is a general result, and can be easily un- 
derstood. The availability of the second charge state provides 
a new non-radiative path for the excited state to decay (via the 
tunneling of an electron onto the partially occupied HOMO). 
This will decrease the probability of the molecule to be in the 
excited state, and therefore decrease the electroluminescence. 
At the same time, the extra charge state also provides what is 
eflfectively an extra conductance channel, thereby increasing 
the current. Both eff'ects will decrease the quantum yield. We 
will come back to this when discussing the results for H2TBPP 
in Sec. IIIB. 
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FIG. 5. (a) Measured conductance of six separate single ZnEtioI 
molecules (A-F) and the oxide surface. [From Qiu et al. (Ref. 51) 
reprinted with permission from AAAS.] (b) Calculated conductance 
of molecule A. The electronic coupling is highly asymmetric (Fs = 
0.637 //eV and Fd = 10 meV) and T = 77 K. 



III. RESULTS 
A. ZnEtioI 

The equilibrium configuration of the ZnEtioI molecule ob- 
tained from DFT calculations is shown in Fig. 1(a), where 
the S4 "saddle shape" structure of the isolated molecule is in 
agreement with the STM topography analysis of Qiu et alJ^ 
Figure 5 shows the measured and calculated conductance of 
the molecule in a STM junction. Qiu et al have measured the 
conductance of six separate single ZnEtioI molecules [A-F in 
Fig. 5(a)], of which only A and B were observed to luminesce. 
These molecules have in common that their conductance plots 
show two peaks within the bias window. This is consistent 
with our rate-equation model, since both the HOMO and the 
LUMO have to be within the bias window for electrolumines- 
cence to be possible, and each will give rise to a peak in the 
conductance. The calculated conductance for molecule A is 
shown in Fig. 5(b). Since only broadening due to temperature 
is taken into account in our rate-equation approach, and not 
due to the coupling to the leads, the calculated conductance 
peaks are much sharper than the measured peaks. It is striking 
that the relative peak heights of the HOMO and the LUMO are 
nearly identical to the measurement, even though the HOMO 
and LUMO are assumed to couple equally to the leads in the 
calculation. The diff'erence in peak height is caused solely by 
the degeneracy of the LUMO and the multiplicity of the ex- 
cited state. 

The measured and calculated bias voltage and current de- 
pendence of the electroluminescence spectrum is shown in 
Fig. 6. The gradual increase in the photo-emission intensity 
at higher bias voltages in the measurement is not clearly re- 
produced in our rate-based calculation, since level broaden- 
ing is not taken into account. Comparison of the calculated 
spectrum to the measurement, on the other hand, shows good 
agreement. Both spectra show several small side peaks be- 
yond the main peak at 790 nm and a series of peaks around 
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FIG. 6. (a) Measured electroluminescence spectrum of ZnEtioI as 
a function of bias voltage. [From Qiu et al (Ref. 51) reprinted with 
permission from AAAS.] (b) Calculated spectrum as a function of 
bias voltage (77 K, assuming instantaneous vibrational relaxation), 
(c) Calculated current dependence of the 790 nm peak intensity at 
2.35 V. Compare to the inset in (a). 



900 nm. The temperature of the experimental setup (77 K) 
makes the identification of individual vibrational modes dif- 
ficult, but the calculation indicates that only a few of the 213 
modes are active in the electroluminescence (see Fig. 7). The 
peak at 790 nm is dominated by the vibrational ground-state 
to ground-state transition, while the peaks at 900 nm consist 
of a handful of modes involving pyrrole breathing and twist- 
ing modes and rotations of the methyl/ethyl side groups. The 
most important modes can be seen in Figs. 7(b)-7(d). The rea- 
son for these modes to be active in the electroluminescence 
spectrum becomes clear when looking at the change in the 
Coulomb potential due to the emission of a photon [Fig. 7(a)]. 
This change is defined as the change in the electrostatic poten- 
tial felt by the nuclei due the diff'erence in the electron charge 
density between the ground state and the excited state. The 
gradient of the potential, and therefore the force, is largest 
where the potential suddenly changes sign [between the red 
and blue areas in Fig. 7(a)]. Comparing the potential gradient 
with the most important vibrational modes shows that those 
atoms move which are close to a large potential gradient. 

Compared to the measurement, the peaks at 900 nm are 
lower with respect to the main peak at 790 nm than they are in 
the measurement. This may in part be caused by the electrolu- 
minescence background due to the NiAl substrate around that 
wavelength, or a varying sensitivity of the charge-coupled- 
device camera in the spectral range, but is most likely mainly 
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(d) 1614 cm- 



FIG. 7. (a) Calculated change in the Coulomb potential felt by the 
ZnEtioI nuclei due to the emission of a photon, [(b)-(d)] The three 
most important vibrational modes responsible for the peaks around 
900 nm in Fig. 6(b). The motions of the nuclei can be directly related 
to the potential gradients in (a). 



caused by the limited number of vibrational quanta taken into 
account in the calculation,^^ resulting in the transitions at 
higher wavelengths to be under-represented in the spectrum. 

In a different measurement, Qiu et al. report observ- 
ing equidistant vibrational features with a peak spacing of 
40+2 meV [Fig. 5C in Ref. 51]. They suggest these peaks 
are possibly higher harmonics of the same vibrational mode. 
In the calculation, however, none of the vibrational modes 
has a sufficiently large electron-phonon coupling^^ to produce 
such a ladder. The calculation does show a series of different 
vibrational modes with a non-zero electron-phonon coupling 
spaced approximately 40 meV apart. 

Besides a dependence on the voltage, Qiu et al?^ also find 
a linear dependence of the photon count on the current [see 
the inset in Fig. 6(a)]. This linear dependence is reproduced 
in our calculations [Fig. 6(c)] and can be easily understood. In 
the case of asymmetric couplings, where the source electrode 
(STM tip) is much more weakly coupled than the drain (sub- 
strate), the average occupation of the excited state is nearly 
zero: it will take a long time for an electron to tunnel onto 
the molecule, but once it is there, it will tunnel off' almost im- 
mediately. The coupling of the source electrode (Fs) can be 
varied by changing the vertical position of the STM tip. This 
changes both the current through the molecule and the photo- 
emission rate, as the latter is directly proportional to the av- 
erage occupation of the excited state. Since the ratio between 
the photo-emission rate and the electron transmission rate is 
independent of Fs (~ see Sec. II), this results in a lin- 
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FIG. 8. (a) Measured electroluminescence spectrum of H2TBPP at 
positive and negative bias. [Reprinted figure with permission from 
Dong et al. (Ref. 52) copyright (2004) by the American Physical 
Society.] (b) Calculated electroluminescence spectrum (assuming in- 
stantaneous vibrational relaxation). The electronic coupling is nearly 
symmetric (Fs = 16.4 yueV and Fd = 4.8 //eV) and T = 300 K. 
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FIG. 9. (a) Calculated change in the Coulomb potential felt by the 
H2TBPP nuclei due to electroluminescence from the LUMO. The 
potential change for the LUMO+ 1 is similar, but rotated by 90° . [(b)- 
(d)] The three most important vibrational modes responsible for the 
peaks around 723 nm in Fig. 8(b). 



ear dependence of the photon count on the current. Note that 
linearly changing the current by varying Fs is only possible 
when Fs <^ Fd. 



B. H2TBPP 

In the measurements of Dong et al the molecule is cou- 
pled almost symmetrically to the leads. The equilibrium con- 
figuration of the H2TBPP molecule used in the calculation is 
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FIG. 10. (a) Measured quantum yield of H2TBPP as a function of 
bias voltage. [Reprinted figure with permission from Dong et al. 
(Ref. 52) copyright (2004) by the American Physical Society.] (b) 
Calculated quantum yield, assuming a 100 % detection efficiency of 
the emitted photons. 



shown in Fig. 1(b). The C2v conformation, and, in particu- 
lar, the angle between the porphyrin center and the four side 
groups, is in agreement with the STM measurements of Jung 
et alJ^ 

The measured and calculated electroluminescence spectra 
at positive and negative bias are shown in Fig. 8. Note that 
electroluminescence at both positive and negative bias is only 
expected when the molecule is (nearly) symmetrically cou- 
pled to the leads (see below). The small asymmetry in the 
coupling causes the photon count at negative bias to be ap- 
proximately half of that at positive bias. As in the measure- 
ment of Qiu et al.,^^ the spectrum consists of a main peak 
around 658 nm, corresponding to the vibrational ground state 
to ground-state transition, and then another peak at around 
723 nm. As with the ZnEtioI calculation, the lower-energy 
transitions are somewhat under-represented in the calculated 
spectrum. The temperature of the experimental setup (300 K) 
again makes the identification of individual modes impossible, 
but the calculation predicts only a few active modes around 
the energy of the peak at 723 nm. These correspond to pyr- 
role breathing and twisting modes — explaining why spectra 
of the two molecules are so similar — and to rotations of the 
side groups (see Fig. 9). 

The measured and calculated quantum yield as a function 
of the bias voltage are shown in Fig. 10. The diflTerence in 
absolute magnitude between the measurement and the calcu- 
lation is due to the unknown detection efi&ciency in the mea- 
surement, which is taken to be a 100 % in the calculation. In 
both the measurement and the calculation, the quantum yield 
is zero until an excited state becomes available in the bias volt- 
age window. It then increases until a maximum is reached 
when all vibrationally excited states have become available. 
As noted above, the quantum yield depends on Fd, and not 
on Fs. Since going from positive to negative bias eflTectively 
means switching Fd and Fs, the maximum value of the quan- 
tum yield changes by a factor of ^ , and is therefore directly 
proportional to the asymmetry in the electronic coupling. 

In the measurement, the quantum yield drops when the bias 
exceeds 3.5 V, which is attributed to damage to the molecules 
in Ref. 52. This can, however, also be explained by the ap- 
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FIG. 11. Calculated current dependence of H2TBPP at positive and 
negative bias. The current is varied by changing the electronic cou- 
pling to the STM tip. 
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FIG. 12. High-energy part of the calculated electroluminescence 
spectrum of H2TBPP at 2.5 V as a function of the vibrational relax- 
ation time, r = 15 K for clarity, the other parameters are the same as 
in Fig. 8(b). Note the appearance of peaks at around 600 nm, to the 
left of the main peak, when the vibrational relaxation time becomes 
comparable to the electron transmission rate. 



pearance of another charge state in the bias window, which 
changes the quantum yield by a factor of ^ 2rs^+rs ^^^^ 
With approximately equal couplings to the source and drain 
(Fs ~ Fd), the quantum yield is reduced by about 50% [as is 
the case in Fig. 10(b)]. However, with very asymmetric, cou- 
plings (Fs ^ Fd), as is the case in the measurement of Qiu et 
al.,^^ the reduction is expected to be unobservable. 

Dong et al.^^ do not report on the current dependence of the 
photon count. However, in the case of symmetric coupling, 
a linear dependence of the photon count on the current, by 
varying Fs with the STM tip, is still expected at positive bias 
(see Fig. 11). At negative bias, on the other hand, varying 
the position of the STM tip equals changing Fd, resulting in a 
decrease in the photon count with increasing current, albeit at 
a lower rate. This is a direct result of the dependence of the 
quantum yield on Fd via 2 ^ ^ (see Sec. II). 

Since vibrational relaxation is several orders of magnitude 
faster than photo-emission (which is on the order of 10"^ s 
in the measurements and calculations under discussion), it 
is assumed in standard fluorescence measurements that the 
molecule is always in the vibrational ground state before emit- 
ting a photon, although at finite temperatures the lowest vi- 
brational excitations can still be occupied, e.g., as observed 
in anti-Stokes scattering. In transport measurements, how- 
ever, even the higher vibrationally excited states are continu- 
ously repopulated by new electrons from the leads. When the 
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electron transmission rate becomes on the order of the photo- 
emission rate, vibrational relaxation will not always take place 
before the electron leaves the excited state. It is therefore the 
ratio between the vibrational relaxation rate and the electron 
transmission rate, not just the photo-emission rate, that deter- 
mines the emission spectrum. This can be seen in Fig. 12, 
where we have plotted the high energy part of the emission 
spectrum of H2TBPP for different values of the vibrational re- 
laxation time. Peaks start appearing to the left of the main 
peak, i.e., at higher energies, as soon as the vibrational relax- 
ation rate becomes comparable to the electron transmission 
rate (at around 10"^^ s). This effect may be present in the 
electroluminescence spectrum measured by Dong et al.^^ [see 
Fig. 8(a)], where such a peak can be seen at around 40 nm to 
the left of the main peak. 



spectra for both symmetric and asymmetric couplings and 
provide an explanation for the suppression of the quantum 
yield at high bias. Although in the measurement of Dong et 
al.^^ this suppression may have been caused by damage to 
molecules,^^ a significant reduction in the quantum yield is 
expected in general for every (nearly) synmietrically coupled 
single-molecule junction in the sequential tunneling regime. 
In addition, we have shown that vibrational relaxation rates 
become important when they are comparable to the electron 
transmission rate, giving rise to peaks in the spectra at higher 
energies higher than the HOMO-LUMO gap. However, a de- 
tailed study of this effect requires a more sophisticated model 
for vibrational relaxation, since a single relaxation time for 
all vibrational modes is only a sufficiently accurate approxi- 
mation when relaxation is either much faster or much slower 
than all other rates. 



IV. CONCLUSIONS 
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